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We propose a scheme to increase the sensitivity and thus the detection volume of nanoscale single 
molecule magnetic resonance imaging. The proposal aims to surpass the Ti limited detection of 
the sensor by taking advantage of a long-lived ancilla nuclear spin to which the sensor is coupled. 

We show how this nuclear spin takes over the role of the sensor spin, keeping the characteristic 
time-scales of detection on the same order but with a longer life-time allowing it to detect a larger 
volume of the sample which is not possible by the sensor alone. 


INTRODUCTION 

Nuclear magnetic resonance (NMR) spectroscopy [Tl|2] 
provides a label-free method for chemical analysis, which 
allows to gain detailed structural analysis as well as dy¬ 
namical information from the specimen under study. In¬ 
deed, e.g. 80% of all structure information on proteins 
stem from the method. However, NMR is rather insen¬ 
sitive and requires chemical purification of the structure 
under study. Much effort has been directed to applying 
NMR to nanoscale samples. Indeed, NMR detection of 
a (4nm)^ voxel of protons has been achieved with mag¬ 
netic resonance force microscopy, a challenging experi¬ 
mental technique operating at ultralow temperature in 
vacuum Under ambient conditions, microcoil de¬ 

tectors have enabled the detection of liquid samples of 
(SOOOnm)^ volume, but further miniaturization of this 
technique is not straightforward 0 0. Magnetic reso¬ 
nance imaging with few or single proton sensitivity would 
have significant impact on a number of research fields in¬ 
cluding life sciences and solid state physics. In addition 
if the detected spins do have good quantum properties 
they might become a valuable resource for quantum in¬ 
formation protocols [?]• 

Over the past decade single nitrogen-vacancy (NV) 
centers in diamond HE] have been proposed and used 
as a novel atomic-size magnetic field sensor for detect¬ 
ing nanoscale nuclear spin ensembles or even single nu¬ 
clear spins [12]. This center is a joint defect in the car¬ 
bon lattice of diamond, consisting of a substitutional ni¬ 
trogen atom and an adjacent vacancy. Its spin triplet 
{S = 1) groundstate can be polarized and read out opti¬ 
cally, so that electron spin resonance experiments can be 
performed on a single spin. A single center can be used as 
a nanoscale magnetic field sensor, able to detect a mag¬ 
netic field in the nanotesla range in an integration time of 
Is [13]. This corresponds to the field of a single nuclear 
spin at a distance of a few nanometers. Indeed, detec¬ 
tion of single nuclear spins, and NMR signals from a 
(5 — nm)^ voxel of various fluid and solid organic samples 
under ambient conditions has already been demonstrated 


[HHE]. In these experiments the detection volume of the 
sample is shown to be limited by the life time of the sen¬ 
sor spin [To] [20] . In this paper we propose a detection 
scheme that is not limited by the Ti of the sensor spin 
allowing for larger detection volumes on a given sample, 
and improving the sensitivity of resolving the distances 
between adjacent spins in the sample. In our prototype 
system, the NV center in diamond, the electronic spin 
interacts with two types of nuclear spins. One type are 
nuclear spins within the sample we want to measure. For 
those spins the electron spin acts as the quantum sensor 
sensing the fluctuations of the magnetic field caused by 
nearby spins in a given sample. Due to the distance de¬ 
pendent (dipolar) coupling between the sensor and these 
spins, the fluctuations of the later are reflected in observ¬ 
ables changes of the sensor spin state. The other type of 
nuclear spins are those in the diamond material itself 
(e.g. or which are much stronger coupled to 

the NV spin, show long coherence times and thus can be 
used as ancilla quantum bits. While the strongly coupled 
spins, close to the sensor inside the diamond have dom¬ 
inant contribution to these changes, the weakly coupled 
sample ones (that are farther away) are mostly masked 
by the stronger ones. To resolve these weakly coupled 
spins in the sample, the strongly interacting spins typ¬ 
ically are dynamically decoupled (DD) [21] so that 
the weak effects caused by farther spins are slowly accu¬ 
mulated in the sensor which is then finally readout m- 
These weakly coupled sample spins can also become un¬ 
detectable when their inverse coupling strength becomes 
comparable to the relaxation time of the sensor. While 
the masking of strongly interacting spins could be re¬ 
moved to a large extent by DD sequences, there is no 
way around to escape the relaxation time of the sensor 
placing an upper bound on the detectable spin volumes 
in the sample. To surpass this problem we propose an 
indirect or inductive method for sensing where a sensor 
spin is coupled to a long-lived ancillary spin (a strongly 
interacting nuclear spin) which now acts as a sensor for 
times longer than Ti allowing one to detect larger parts 
of the sample. 

The physical setup for our proposed scheme consists of 
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FIG. 1. a Schematic representation of the central idea, depicting the increase in detection volume of a sample on top of the 
diamond both through direct and indirect sensing, b The energy-levels of the sensor-ancillary system where only two levels 
of the sensor (within the triplet subspace) are resonantly coupled to a two-level ancillary (nuclear) spin, c We show the pulse 
sequence for a direct detection of a spin by the sensor. In the same hgure we show the sharp transition of a sensor observable, 
{Sx) at T = + J). The parameters chosen for the simulation are ujl!J = 10, ^ = 7r/4 and the optimal number of 

repetitions N — 2Z. d Mapping pulse sequence Um(jq) through which the sensor-sample coupling is directly reflected on the 
dynamics of the ancillary spin F. Also shown in the hgure is the variation of the observables {Sx) and {Fx) both in the presence 
and absence of mapping. 


a sensor described by a three-level (spin-1) system that is 
interacting with the spins (that are typically spin-1/2 )of 
the sample. In the reference frame of the sensor one can 
think of the spins in the sample that are all distributed 
in a three-dimensional plane comprising of the NV axis, 
(the direction along which the sensor is aligned) and the 
plane perpendicular to it. For simplicity we consider the 
two-dimensional conhguration where the position of each 
spin in the sample is described by polar coordinates r 
and 0. We assume that the sensor is placed along the 
2 ;-axis and is parallel to the applied magnetic field. In 
addition to the spins of the sample we also have an an¬ 
cillary spin which is strongly interacting with the sensor 
along the ^-axis. The relevant part of the Hamiltonian 
that determines the dynamics is simply given by 

N 

^ Jf + gS^F^ + ^ Ji (cos eat + sin . 

i i=l 

( 1 ) 


where the sensor, ancillary and sample spins are respec¬ 
tively denoted by S', F and I. The zero-field splitting 
of the sensor and the magnetic field strength Bq can be 
chosen such that the sensor can act as a two-level sys¬ 
tem formed by the magnetic levels = 0, — 1 (see 

Fig. 1(b)). In the presence of field Bq the sample spins 
precess at a Larmor frequency uol- 


INDIRECT SENSING BY AN ANCILLARY SPIN 

In a standard detection scheme, using the sequence 
shown in Fig. Ic, mi the sensor is repetitively flipped 
(dynamically modulated) about the plane perpendicular 
to the z-axis at a rate that matches the energy of a par¬ 
ticular sample spin li which is = 2uol F Ji cos Oi M- 
Under such matching conditions there is a transition in 
the sensor spin state which leads to a large observable ef¬ 
fect in its optical or electrical readout (see Fig. Ic). The 
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peak positions r = depend on the parallel compo¬ 
nent {z — z) of the dipolar coupling between the sensor 
and the sample spins, while the width depends on the 
perpendicular component {z — x) of the interaction i.e., 
A « N/JisinOi, where N is the total number of pulse 
repetitions (see Fig. Ic). The optimal number of repeti¬ 
tions Nopt required to obtain maximum contrast (height 
of the peak value) again depends on the values of Ji and 
Oi. We will show later that in the presence of dissipa¬ 
tion the peak width increases and the height decreases 
(thereby increasing the uncertainty of the detected posi¬ 
tion) but the peak position remains unchanged. Hence, 
by switching on a radio frequency (RF) field along x- 
direction one can exchange the roles of parallel and per¬ 
pendicular components of dipolar coupling so that the 
peak now appears at = 2uul-\- Ji sin Oi. Thus both the 
variables can be readout through the peak-positions only. 
This detection method can also be understood as the well 
known spectral overlap of the applied control (the filter 
function) on the sensor with the spectral distribution of 
the sample spins (Fig. Ic), a widely used technique in 
magnetic resonance imaging. 

As can be seen, the above detection scheme can only 
resolve spins that are spectrally separated by distances 
larger than the width of the detection peak. Over the 
past years there have been quite some advances in im¬ 
proving this resolution by using more complex pulse se¬ 
quences so as to have better resolution by avoiding cross¬ 
talk between the nearby spins, and also a high degree 
decoupling to the unwanted dephasing noise spectrum. 
But, the unavoidable relaxation of the sensor still remains 
a barrier, placing an upper bound on the detectable vol¬ 
umes of the sample. In the current proposal we sur¬ 
pass this effect by using a long-lived ancillary spin that 
is strongly coupled g Ji to the sensor. 

To make this ancillary spin which is not directly cou¬ 
pled to the sample spins also an indirect sensor, the first 
step would be to map the sensor-sample dynamics onto 
the ancillary spin. Though the below discussion is valid 
for any number of spins in the sample, for simplicity 
we shall concentrate on the case where the sample is 
described by a single spin coupled to the sensor with 
strength J. The dynamics of this three-spin system, the 
central spin (sensor) is coupled to the sample spin I and 
to the ancillary spin F is analytically solvable. Before we 
go ahead to design a mapping pulse sequence, we would 
like to first evaluate the behavior of the ancillary un¬ 
der free evolution generated by the Hamiltonian given in 
Eq. (1) as a function of J. One finds that the the time- 
dependent evolution of the coherence for the ancillary 
spin, given by 

Cf{t) = (Fa;) +i(Fj^) = i[l + cos(flT)]CF(0), (2) 

has no dependence on the sample parameters J and uj 
for any time r. This can be seen from the structure of 




FIG. 2. Using the detection sequences Ds and Dp shown on 
the top of the figure we plot the expectation values of the 
sensor and ancillary spins {Sx) and {Fx) as a function of the 
Larmor frequency ul both in the presence and absence of re¬ 
laxation. The sensor-ancillary coupling is much stronger than 
the sensor sample spin coupling and is taken to be ^ = 80 J, 
and 0 = 7r/4. Relaxation is introduced through the Krauss 
form given in Eq. (5) with jsjg — 10“^. For consistency we 
have considered the relaxation of the sensor during the ancil¬ 
lary control pulses in the Df sequence as they are assumed to 
be considerably slower than the manipulations on the sensor 
spin. 


the Hamiltonian (Eq. 1), where the ancillary-sensor cou¬ 
pling term commutes with the remainder and hence the 
sensor-sample coupling {S — I) only leads to a overall 
phase shift and does not contribute to the dynamics of 
the ancillary spin. The situation becomes different if a 
microwave field that couples energy states of the sensor 
is switched on. The total dynamics cannot be separated 
into the dynamics on S — F and S — I subspaces as these 
terms no longer commute resulting in direct correlations 
between the ancillary and the sample. Instead of switch¬ 
ing a continuous microwave field, we rather do it strobo- 
scopically using the pulse sequence shown in Eig. 1(d). 
Under the action of this pulse sequence, the coherence 
of the ancillary spin now oscillates as a function of J as 
shown in Eig. 1(d). Eor a fixed time tq = it! g and Hq = 0 
the analytical expression for Cp takes a simple form 


Cf{t = To) = cos(Jro) — i cos(6>) cos(Jro/2) sin(Jro). 

(3) 

Now one can clearly see that on time-scales of tq, the 
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FIG. 3. Using the detection sequence Dsf shown on the top 
of the figure we plot the expectation values of the sensor and 
ancillary spins {Sx) and (iU) as a function of the Larmor 
frequency lul both in the presence and absence of relaxation. 
The sensor-ancillary coupling is much stronger and is taken to 
he g = 80 J, and 0 = 7r/4. Relaxation is introduced through 
the Krauss form given in Eq. (5) with = 10“^. 

coherence of the ancillary spin purely oscillates on the 
frequency of the sample coupling parameters J and 0. 
For J g the purity of the ancillary state is close to one 
during the evolution, and changes as if it were directly 
coupled to the sample spin /. Equally the coherence of 
the sensor spin on this time scale only depends on J, 
given by 

Cs{ro) = {S,)^i{Sy) 

= ^(1 + cos(Jro)^) — icos{0) cos(Jro/ 2 )^ sin(Jro). 

( 4 ) 

We would like to note that by using the above mapping 
sequence the time units are now renormalized in multi¬ 
ples of To- For example, with a Larmor frequency ujl and 
interaction J, a detection sequence shown in Fig. 1(c) 
would give a sharp transition in the sensor coherence at 
T = 7rl2ujL J cosO. But with the mapping sequence, r 
should now be varied in steps of tq, implying that spins 
whose total Zeeman energy is smaller than the S — F cou¬ 
pling strength g are detectable and the parameters should 
satisfy g > ujl > a commonly available situation in 
real physical systems. Having mapped the parameters of 
the sample onto the ancillary spin, the next step is to 
detect the sample spin by the ancillary. One should also 
remember that the sample spin could also be detected by 
the sensor itself as its coherence is also a function of J and 
0 during every mapping sequence. One can see from Fig. 

2 that we now have the possibility to sense the sample 


spin either by running the detection sequence on the sen¬ 
sor Ds or a detection sequence on the ancillary spin 
and see an almost identical transitions in their observ¬ 
ables {Sx) and {Fx) respectively. Though the difference is 
quite small one can find that sensing on the ancillary has 
a better contrast than on the sensor, and this amplifies 
in the presence of relaxation which we shall discuss later. 
The advantage of the mapping sequence becomes more 
evident when the sensor needs to be reinitialized quite 
frequently, making it ineligible to sense the weak sample 
spin, but the ancillary spin can still reflect the same tran¬ 
sition independent of sensor’s reinitialization rate in the 
absence of any relaxation. If we now squeeze both the 
mapping and detection into one sequence, allowing for 
control only on the sensor (as they could be quite faster 
than controlling the ancillary) but still achieving the goal 
of reflecting the information of the sample spin on both 
the sensor and ancillary could be quite useful. For this 
we use a single asymmetric sequence (shown in Fig. 3) 
where a 7r/2-pulse (for mapping) and tt (for detecting) 
are performed on the sensor during the first and second 
half of the sequence. Surprisingly, repeating this single 
sequence leads to identical dynamics on both the sensor 
and ancillary in the absence of relaxation. This can be 
seen from Fig. 3 where we have plotted the observables 
(Sx) and (Fx) that overlap perfectly. 

We now turn to a realistic situation where the sensor 
is relaxing to its thermal equilibrium at a rate 75 . As¬ 
suming a pure Markovian decay, we introduce relaxation 
as a depolarizing channel [23] through the Krauss form 
such that the sensor spin state at any time t is simply 
given by 

^ (1 - 1 0 Tvsipit)) + ( 5 ) 

where p{t) is the density matrix of the total system at any 
time t. Though the above way of introducing relaxation 
into the problem is simpler than solving the master equa¬ 
tion for the total system exactly, it still could capture 
most of the physics relevant to the current discussion. 

From the above analysis we have already seen that in¬ 
direct sensing through the ancillary spin shows better 
contrast in the presence of relaxation even though they 
have equal contrast in the absence of relaxation. We 
now want to know if the detection through ancillary spin 
could be made independent of the sensor relaxation 75 . 
We have already seen from Fig. 2 that using the map¬ 
ping sequence, an independent detection of the sample 
spin can be made through the ancillary using Dp. While 
running the detection sequence on the ancillary even if 
the sensor is reinitialized though it destroys the sensor- 
ancillary {S—F) correlations, the ancillary-sample (F—I) 
remain intact. Hence the availability of a fully coher¬ 
ent sensor spin state at the end of each repetition would 
map the sample spin parameters onto the ancillary with 
almost unit efficiency (as in the ideal case of 75 = 0 ). 
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FIG. 4. The maximum contrast of the ancillary spin observ¬ 
able is plotted as function of the relaxation time for different 
reinitialization rates for two cases (a) when the sensor reini¬ 
tialization does not dephase the ancillary spin and (b) when 
the sensor reinitialization dephases the ancillary spin at a rate 
Tn- The dephasing Tat is chosen such that = 10“^, 

where the sensor-ancillary coupling and other parameters are 
similar to that given in earlier plots. 


Instead if the reinitializations are performed at a slower 
rate, the relaxation effects of the sensor enter the dynam¬ 
ics of the ancillary, thereby reducing its contrast from the 
ideal value. This is shown in Fig. 4(a), where we plot 
the contrast observed on the ancillary spin readout as a 
function of the relaxation time Ti = l/ys, for different 
initializations of the sensor. Under very frequent initial¬ 
izations the contrast remains independent of Ti, but for 
slower rates the loss of purity occurred till the reinitial¬ 
ization decreases the overall contrast. In a realistic sit¬ 
uation one should also consider an additional dephasing 
incurred by the ancillary during the reinitialization of the 
sensor. This happens for example due to the difference 
in the hyperfine coupling between the sensor and the an¬ 
cillary spin in the ground and excited states. Taking into 
account this additional dephasing, one immediately finds 
that frequent reinitialization may not be the ideal situa¬ 
tion, which can be seen from Fig. 4(b). 

We now make some numerical estimates on the new 
range of distant spins that the ancillary spin could sense. 


For NV centers in diamond one can find an ancillary 
spin such that the sensor-ancillary coupling is ^ ^ IMHz, 
and the sensor has a typical relaxation rate 75 = IkHz. 
The sensitivity with which a sample spin whose coupling 
strength J = lOkHz could be detected reduces to 17% 
of its maximum achievable case when 75 = 0. Now run¬ 
ning the detection sequence on the ancillary spin, with 
the possibility to make a tt flip at rate we find that the 
sample spin could now be detected with an increased sen¬ 
sitivity of 25%. Now if we reinitialize the senor spin at a 
rate of ^/lO i.e., the sensor is reinitialized three times in 
the entire sequence at very 10 /is, increases the sensitivity 
to 33%. The dephasing incurred by the ancillary (nu¬ 
clear) spin during the reinitialization process is propor¬ 
tional to the reinitialization time of the sensor (electron) 
spin Tr and the hyperfine coupling between the ancillary 
and the sensor in the excited state Qq. While choosing 
strongly interacting nuclei as ancillary spin will improve 
the spectral resolution of the sample spins, it equally will 
get dephased fast enough during frequent reinitialization 
of the sensor and vice versa. Thus finding an optimally 
coupled ancillary spin would become a key requirement 
for the indirect sensing proposed in this paper. 


CONCLUSION 

In conclusion we have proposed a new method to detect 
distant spins that are very weakly coupled to the sensor. 
When this coupling becomes much smaller than the re¬ 
laxation rate of the sensor, direct sensing yields almost no 
contrast, while an indirect sensing through a long-lived 
ancillary spin coupled strongly to the sensor can be used 
to detect the same spins with a high contrast. We have 
shown that this kind of dual sensing could be achieved 
by a mapping sequence between the sensor and the an¬ 
cillary so that the information related to the sample spin 
is stored in both the spins, and becomes detectable by 
either of them. Though the analysis is made on a simple 
model system having a single sample spin, the analysis 
could be straightforwardly generalized to highly dense 
samples with few hundreds of the sample spins. 
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